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Energie libere in viscoelasticità e applicazioni 
alle equazioni alle derivate parziali.* 


Mauro Fabrizio! 


1 Introduzione 


Nei materiali con fading memory, lo stato è dato dalla storia. In questo lavoro 
mostreremo come sia più naturale definire lo stato mediante la risposta for- 
nita dallo stress quando è sottoposto ad un processo di deformazione costante. 
Inoltre metteremo in evidenza il significato e l’importanza di determinare una 
norma naturale nello spazio dei processi e degli stati per un materiale con fading 
memory. A tal fine mostreremo come l’idea considerata da Gentili in [12], di 
introdurre una topologia nello spazio dei processi attraverso la continuità del 
lavoro, fornisca una topologia naturale nello spazio degli stati, che Gentili di- 
mostra essere duale rispetto a quella dei processi. Questo punto di vista sembra 
particolarmente utile nello studio dei PDEs connesse a tali problemi. Per questa 
ragione l’articolo finisce con un’applicazione alle equazioni integro-differenziali 
della viscoelasticità. Per tale problema, grazie a questi nuovi spazi e senza al- 
cuna ipotesi sul valore del coefficiente G'(0) e sulla regolarità del nucleo G”(s), 
possiamo provare l’esistenza, l’unicità e il decadimento asintotico delle soluzioni. 


2 Fading memory e termodinamica 


Un materiale viscoelastico è definito da una equazione costitutiva che lega il 
tensore dello stress T' al gradiente di deformazione F mediante un funzionale 
del tipo - 

T(2,t) = Î(F‘(2)), 


dove F*(x,s) = F(x,t—s)èla storia del gradiente di deformazione F. Nel caso 
lineare abbiamo 


T(x,t)= Go(x)E(r,t)+ f G'(x,s)E*(x,s) ds, 


*Questa ricerca è stata realizzata con il contributo del GNFM (INDAM))) e in parte sostenuta 
dal MIUR (Cofin 2002) attraverso il progetto di ricerca “Mathematical models for materials 
science” e dall'Università di Bologna - Fondi per progetti pluriennali. 
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Vu+ (Vu)? 


dove E = ———_ è il tensore di deformazione. Inoltre, Go(x) e —G'(2, s) 
sono tensori simmetrici e definiti positivi, così 


Goo(1) := Go(x) + 1° G'(x,s)ds. 


Naturalmente, questo studio può essere applicato all’equazione costitutiva non 
lineare del tipo 


(e.°) 
T(x,t) = Go(x)®(F(2,t)) +f G'(x,s)®(F*(2,5)) ds, (1) 
0 
dove ® : IR? + IR3 è una funzione dispari rispetto ad F, non lineare e regolare. 


Una applicazione P : [0,d,) — Lin continua a tratti e definita come: 
P(t) “= É(t), LE (0, dp) 


è chiamata processo cinetico di durata d, € IR*. 


ES 


Definizione 1 Un materiale semplice è dato dall’insieme (I1,L, p,T): 


a - Il è l’insieme dei processi cinetici definito come: Il = {P : [0,d,) —> Lin 
continui a tratti; se P,, P. € II, allora P,* Pa € II, dove 


P,(t) se t€0,d,,) 

P, * P.(t) = ; ToARi ; 
1 + Pa(0) { Pa(dp, —t), se t € [dp dp, + dpa) 
inoltre, se P € II, allora Pt, t,) E II, [t1,t2) C [0,dp)} 

b- Xè lo spazio degli stati, i suoi elementi sono indicati con 0, 


c - la applicazione dp: XL xII — X è chiamata funzione di evoluzione, tale che 
se c' è lo stato inizale e P è un processo: p (c',P)= of, 


d- la applicazione T:LxI + T è la funzione risposta che, ad ogni stato 
o e ad ogni processo P, assegna il processo del tensore stress TP su tutto 
l’intervallo temporale [0, d,). Cioè 


TP(t)=T(c,P)(t), tc 0,4). 


Definizione 2 Il sistema è chiamato causale, se la priprietà d della Definizione 
1 si può sostituire con la nuova condizione 


d’ - la applicazione T': X x Lin —> Sym è tale che 


T(t) = T(0(6), P(6)). 
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Osservazione 1 In questo articolo consideriamo solo sistemi causali, per i 
quali, poichè o(t) = p(0, Pi), la funzione T è legata a T dalla relazione 


T(p(0, P.), P(t)) = T(0, P)(t) 


Un materiale viscoelastico è un materiale semplice dove lo stato è dato dalla 
storia o = (E*) = (E(t), Et)e il processo da 


P(t) = É(t), t € [0,d,). 
In seguito consideriamo il lavoro W del processo P, definito come 
do . 
W (0, P) =] T(o(t), P(t)) - L(t)dt. 
l) 


Definizione 3 (Noll, 72) Due strain storie Et, Ei si dicono equivalenti se 
per ogni processo EP : [0,dp) — Sym, soddisfano 


î(E,E?)=f (2,6), 


o analogamente, come provato da Gentili, [12], due storie si dicono equivalenti 
se per ogni processo EP : [0,dp) —> Sym abbiamo 


W (21,2) =W (EE?) . 


Teorema 1 (Del Piero, Deseri, 96, 97) Nel caso lineare, due storie Ei, Ei 
sono equivalenti se e solo se Ex(t) = Ea(t) e 


Sf G'(s + 7)Ei(s)ds = f, G'(s+7)E5(s)ds, Vr>0. (2) 
0 0) 


Secondo la definizione di stato, due coppie (Ei(t), E{(7)) e (E2(t), Ei(7)) 
equivalenti (nel senso di Noll) sono rappresentate dallo stesso stato 0(t) (vedere 
[14]). In altre parole lo stato può essere rappresentato dalla coppia 


o(t)=(E(t),Ît(7))=I*(7), 
dove 


Î(r7)=- ul G'(s +7)E'(s)ds. 


Definizione 4 (Gentili, 02) Un processo ÈP : [0, 00) > Sym è detto processo 
a lavoro finito se 


dp ® . 
W(0', EP) = J T(0, EP) ÈP(7)dr < co. 
0 
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Inoltre, il lavoro W(0î, ÈP ) nel caso lineare si può scrivere nella forma 
wo, EP) = 36wE(dp)- Edp) 
4/ fo G(|r — s|)EP(r) - EP (s)drds (3) 


= 30 oE(dp) P(dp) +3 -f Go(w)EP (w) - EP (w)dw, 


dove G = G(8) — Goo, Ge e ET sono le trasformate di Fourier di G e di ÈP. 
L’insieme dei processi a lavoro finito è dato da 


Ho = {er Pag f° G(|r — s|\.EP (7) - EP (s)drds < co) © 


Questo insieme è uno spazio di Hilbert se il kernel G soddisfa la Seconda Legge 
della Termodinamica (i.e. G.(w) > 0 ) e la norma è data da 


e] 


f È G(|r — s|) EP (7) - EP (s)drds 
o Jo 


z ni n Go(w)EP(w)- EP(w)dw. 


Il dominio di definizione degli stati ammissibili è l’insieme di tutte le storie 
di deformazione che rendono il lavoro ben definito quando il processo appartiene 


a Ho. 


W(I', EP) n 3GE(dp) - E(dp) + 5 T i G(|IT - s|)EP(7) . È? (s)drds 


+G21*(0) - E(dp) — T Î*(r) -EP(s)dr < 00. 
(5). 


Quindi l’insieme degli stati ammissibili Î*(-, E*) è dato dal duale HG di Ho, 
cioè 


Ha = {raf Î(r)- ÈEP(s)dr < co, VEP e Ho} 
{rx if È e*(Ir- 7) - È (7)drdr <co), 


dove (G*)p = G71(w). 


L’insieme Hi è uno spazio di Hilbert se la norma è data da 


fl 


Sf È G*(|r - 7'|)Î*(7) - Î*(7')drdr' 
o So 
1/ cri) Td 
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Definizione 5 Una coppia (o, P) € X x II è detta ciclo se: 
p(o,P)=o. 


Seconda Legge per i processi isotermici. 
Per ogni ciclo(o, P) € LL x II vale la disuguaglianza 


d 
f TemPm). 10 dre>ò (6 


Per i materiali fading memory, i cicli sono molto rari. Per questa ragione con- 
sideriamo una nuova forma della Seconda Legge. 


Forma Forte della Seconda Legge (per i processi isotermici) [7]. 
Sia Ly := {o;3 Pe Il; d(o,P)=0"}. L'insieme dei lavori che si ottiene 
passando da uno stato o ad un arbitrario stato o € Lo 
W (0) := {W (o, P); Pe I} 
è limitato inferiormente. C’è uno stato co , chiamato stato zero, tale che 


inf W (00) = 0. 


Questa ultima definizione della Seconda Legge soddisfa il Principio di Dissi- 
pazione dato da Gurtin & Herrera in [15]. 
Come mostrato in [10] dalla Seconda Legge della Termodinamica abbiamo 


Gi(w) = I° G'(s)sinwsds<0, per ogni w € Rt+ (7) 
mentre dalla Forma Forte della Seconda Legge, si ha 
Goo > 0. (8) 
Definizione 6 Una funzione w : Dy + R* è chiamata energia libera se 


a - il dominio Dy CD=Y è tale che per ogni o € Dy e PE II, lo stato 
o=p(0,,P)€ Dy, 


b- co € Dy eyw(c0)=0, 


c - per ogni coppia 0), 07 € Dy e PE II tale che p(0,,P)= 072: 


Y(02) - 4 (01) < W (01, P) (9) 


Nella viscoelasticità lineare, ci sono molte energie libere [2], [4], [11]. La 
famiglia 7 delle energie libere è un insieme convesso. Questo insieme 7 presenta 
un minimo e un massimo %,n, Vm- 


L'energia libera massima è stata considerata in [7] 
1 
VW (E') = 3GoE (1) E (1) 
1 (o.°) (eo) 
43/, i G12(81 — 82) (E° (81) — E (t)) - (E* (s2) - E (t)) ds2ds9 (10) 
lo) 


0? 
dove G12(s1 — 82) = Deda —— G(|s1 — s2|). 


Se G' (8) <0eG”(s) >0, c'è un'energia libera intermedia, chiamata energia 
libera di Graffi-Volterra, data da 


We (E') = 50xE (t)-E(1) 
-;/c0E-50)-(E0-£0) 4 1) 


L'energia libera minima è stata trovata da Bruer e Onat nel 1964 [1] nella 
seguente forma 


Vm(ÉEm ) = 50E(t)- -E(t +3/ Le G (|s1 “= 52) Em(s1)- Ban (52) ds1ds2, 


(12) 
dove È, è il processo ottimale che fornisce il massimo lavoro recuperabile, ma 
non è un funzionale della storia E* (s) o dello stato I* = (E(t), I*(7)). 

Golden [13] ha dato una rappresentazione dell’energia libera minima in ter- 
mini di /*(r) come 


Um(') = 30% E(t)-E(t)+= fo 1 H(r,1)ft(r)-Ît(r)ardr!, (13) 


dove H(7r,r') è un opportuno kernel dipendente da G. 3 i 
Introduco infine alcune nuove energie libere scritte come funzioni di /*(.). 
La prima è data da 


Vi(#) = 30oE@-E0+3 3/_ 1a G*(ir = )E (7) Er')drdr' 
= 3OwE@)-EM+3/" GT) 0)da 
dove (G*)r = G-! > 0. La seconda è 


v(Î) - 3GE(0) - E(t) + 5f G-1(7)Î*(7) - Î*(r)dr 


dove G(s) < 0, G(s) > 0. 
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3 Applicazioni alla viscoelasticità lineare 


3.1 Comportamento asintotico 


Consideriamo l’equazione alle derivate parziali connessa con un materiale vis- 
coelastico e definita nel dominio Q= Q x RT 


i(x,t) =V-(Go(x)Vu(2,t) + È G'(x,5)Vu'(x, s)ds) 
— (4) 
=V-(Go(x)Vu(2,t) +/ G'(r,s)Vu'(x, s)ds) +V-Î°%x.t), 
dove Î°(,t) = /s° G'(x,8 + t)Vut=%(x, s)ds. Assumiamo inoltre come con- 
idizioni iniziali 
u(7,0) = uo(2), (7,0) = o(x) ; (7,5) = (7,5), se Rt (15) 


e come condizioni al contorno 
Uan = 0; (16) 


Per ottenere una definizione rigorosa della soluzione debole nell’intervallo 
temporale IR*, dobbiamo introdurre lo spazio di funzioni 


G(R*, V(0) = {u € Lie(IR, V(09)); 
f ta (G(7 — T)u(7), u(7'))y drdr' < co} 
e il suo duale G'(IR*,V(9)). Indichiamo poi 


£(IR*, H3( ci {u e G(IR*, Hi(9)); è e GAIR+, L?(Q )} 
L*(R+, H3(0)) = {ue G(IR+, AL(O)); ù e G(R+t,L2(0))} 


Ora consideriamo il nuovo spazio £(IR, H}(Q )) ottenuto dall’insieme delle trasfor- 
mate di Fourier pe per ogni u € £(IR*, H}(Q)) e indichiamo con l’(IR, H}(9)) 
il duale di É(IR, H3(9)). Sicuramente il teorema di Plancherel per le trasformate 
di Fourier definisce un isomorfismo naturale tra £(IR+, H}(9)) e É(IR, Hi(O)). 


Definizione 7 Una funzione u € L(IR*, Hi(Q)) è detta soluzione debole del 
problema differenziale (14)-(16) con uo € H*(Q), do € L?(9), e u%(x,-) tale che 
Ît=0(£,7) = = fa” G'(x,5+7)Vu(7, s)ds € L'(IR+,L?(0)), se u(2,0) = uo(2) 


quasi ovunque în { e 
“fa } _ u(r ‘@' s)Vu! . ©. 
i VEGRRCO) {cov ( + f G'(s)Vu'(2, s)de} Vép(x.t)dx 
= -/ to(x) - $(7,0)dr + J ÎI*=%(x,7)- V@(x,r)drdx (17) 
Q 0 a 


per ogni $ € L*(IR*, Hi(9)). 


Se indichiamo con a(u, $) la forma sesquilineare su £(IR+, Hi(0)) 
‘00 t 
a(u, dg) = J / u(c,t)-d(2, )-{Gxvue, t) + G'(s)Vut(x, s)ds}-VA(c.d)de 
0 Q 0 


per ogni $ € L*(IR*, Hi(0) e tale che 9(,0) = 0, allora l’equazione (17) si 
può scrivere come 


a(u, $) = -J to(x) - $(,0)dr +/ / Ît=%(x,7)- V(x, r)drdr. 
Q 0 Q 
Adesso possiamo enunciare il seguente 


Teorema 2 Sotto le ipotesi (7)-(8) per la funzione rilassamento G, il problema 
(14)-(16) ammetté una ed una sola soluzione debole u € L(IR*, Hi(OQ)). 


Con un semplice cambiamento di variabili, è sempre possibile ottenere la con- 
dizione iniziale uguale a zero. Quindi, senza perdere in generalità supporremo 
uo = 0, to = 0. Indichiamo con è la forma sesquilineare su £(IR, H3(0)) 


00 
a(d,0) = - J , utt ion 
— 00 


+2 f n / [CO (2) +6 (2,0)] v'è (210) ve' (2,0) ded 


per ogni $ € É(IR, H}(9)). Il teorema di Plancherel applicato a (17) dà 
al@= = fo fit) we (7,10) ded (18) 
27 -—0 SL 


Lemma 3 Una funzione è e É(IR, Hi(Q)) è una trasformata di Fourier di 
una soluzione debole del problema a valori iniziali a contorno (14)-(16) nel 
senso della Definizione 5 se e solo se vale l’uguaglianza (18) per ogni d E 
L'(IR, HI(9)). 

Ponendo nella (18), d (7,7) = 41 (1) 49 (w), con p; € H!(I) e p, € L?(R), 
dalla scelta arbitraria di 4, segue che, per quasi tutti w € IR, vale la seguente 
identità: 


fa (ass) Gi (a+ | [GI (2) +6 (21 u)] vd (esi) wRi (e) de 
a Q 
- fit) vale, (19) 


per ogni 9, € H*(Q). Ma l’identità (19) significa che è (-,w) è una soluzione 
generalizzata in H! (Q) per il problema ellittico 


ud (2,1) -7-{[60(2) +0 (z,w)] vd (2,6)} = v-Ît=% 2,0) 
ù(cr,w)=0, rEd0. 
(20) 
Seguendo la dimostrazione del Teorema 3 di [9] possiamo provare 
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Lemma 4 Per ogni w € KR il problema (20) ammette una e una sola soluzione 


ù(-,5w) € H!(Q). 
Dalla parte reale di (20) abbiamo 
w |a (c,w) <K Iva (c,w)+|y TE 0) è (2,0), = (21) 


dalla parte immaginaria di (20) otteniamo 


/ Ge (2,0) lov è (2,0) de = J WE, (0) leg è (2,0) de 
e) a 


pi 


Infine, sotto le ipotesi sulle condizioni iniziali, la disuguaglianza (22) e l'equazione 
(201) in un intorno di w = 0 forniscono la disuguaglianza 


|a + va], <e f, 


Per cui da (23) e (21) abbiamo: 


questo comporta che è € É(IR, H}(9)) . Inoltre l’isomorfismo tra £(IR+, Hi (0) 
e L(IR, H}(9)) garantisce che è sia la trasformata di Fuorier della soluzione 
u € L(IR*, Hi(Q)) del problema (14)-(16) . 


IA 


Ît°(z,0)-wvî (x, )| dx. (22) 


I° (2,0) wyd (z,0)| dr < 00, (23) 


|, < 00, 


CÈ (w)(1+ |o|)Wa (w)|l + 


GI (0)(1+ lo] (0)[[, < © 


2 
L 


3.2 Teoria dei semi-gruppi 


Nello studio del sistema differenziale (14) con la teoria dei semi-gruppi, possiamo 
scrivere il problema (14) nella forma 


u(c,t) = v(2,t) (24) 
Ù(£,1) V.(CoVu(2,8)+ / G'(x,s)Vui(x, s)ds) (25) 


Ei(z,s) = -L.Ei(,9)- É(2,1) (26) 


dove Et(x,s) = E'(x,s)- E(x,t). Oralo stato è dato da x = (u,v, Ej) € K. Nei 
lavori [3], [8], [9] lo spazio X è un sottoinsime del dominio Hj(Q x L*(Q)x Sy, e 
il teorema sulla stabilità è stato dimostrato con le condizioni iniziali nello spazio 
K. 

Quando usiamo la funzione Î*(r), il sistema (24)-(26) è equivalente al prob- 
lema 


ù(a,t) = v(£,1) 
o(c,t) = V-(GwoVu(z,t)+1Î*(0) 
Ît(x,7) = Lie) - (MPA) 
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dove lo stato x = (u,v,Î '). Adesso possiamo usare l’energia libera yy come 
una norma sullo spazio degli stati, in questo modo è facile provare che lo spazio 
degli stati coincide con il dominio Hj(2 x L?(9) x Sy,, dove Sy, D Syo: 


Teorema 5 Sotto le ipotesi 
1. C'(w)= So G'(s)sinwsds < 0, w>0 
2. Ja e Rt+;0<-G"(s) < aG"(s), s e Rt 
ogni soluzione (u,v, Î*), con le condizioni iniziali xo € Sy,, è tale che 
(Ilù(E)IT= + Pa(Vu(t))) < Me *(Ilt(O)]{72 + V2(Vu(0)) 


dove M e u sono delle opportune costanti. 
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